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Abstract 

We propose a new method for constraining a primordial homogeneous mag- 
netic field with the cosmic microwave background. Such a field will induce an 
observable parity odd cross correlation between the polarization anisotropy 
and the temperature anisotropy by Faraday rotation. We analyze the neces- 
sary experimental features to match, and improve, current constraints of such 
a field by measuring this correlation. 
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In our Galaxy, as well as other spirals, large scale magnetic fields on the order of 10~^ 
Gauss have been observed, the origin of which is unknown Schemes that invoke a dynamo 
mechanism to explain these fields all rely on the presence of an initial seed field [Q], prompting 
numerous explanations ranging from nonlinear battery mechanisms, to changes in the nature 
of the electroweak force 0. Furthermore, the mechanism of dynamo generation itself has 
yet to be developed in a fully self-consistent manner, and the maximum possible fluxes may 
have been overestimated Related to the origin of magnetic fields in spiral galaxies is 
the question of the origin of the 10~^G fields detected in high redshift galaxies and in 
the damped Lya clouds . At these early times a plausible dynamo explanation has yet to 
be proposed. Alternatively, large scale galactic and extragalatic fields can be explained by 
the adiabatic compression of a primordial magnetic field on the order of 10~^ Gauss today 
0. Thus, while it is certainly possible that several of the physical schemes being studied 
may play a role in the formation of large scale magnetic fields, the considerable debate 
surrounding this subject can not be resolved until the initial magnetic field configuration is 
better known. This can only come from further observation. 

The cosmic microwave background (CMB) supplies us with the oldest and most extensive 
probe of the early universe. It has recently been argued that it may provide tight constraints 
on primordial magnetic fields. In |Q] it was shown that a random magnetic field with an 
amplitude of order 10~^G would lead to distortions in the angular power spectrum of the 
CMB of around 10%, observable with planned satellite experiments. In |^ it was shown 
that, if one included the full anisotropic effects of the magnetic field one could use the 
existing large angle measurements of the COBE satellite to set a limit on a homogeneous 
primordial magnetic field today, Bq < 6.8 x 10~^(fio^^)^^^ G where f2o < 1 is the cosmological 
density parameter and the Hubble constant is Hq = lOOh km s^^ Mpc~^. Other authors 
have suggested that Faraday rotation might affect the CMB anisotropy and polarization in 
a distinct way. In |jlO|, the authors derived analytic expressions for its effect on the dipole, 
quadropole, and octopole of the CMB anisotropy in a Bianchi I cosmology. More recently 
in illj, the authors found the effect that a homogeneous field with an amplitude of 10~®G 



would have in decreasing the polarization of the CMB, while in |12[ it was proposed that a 
two frequency measurement of the polarization should be able to set bounds on a random 
magnetic field on the order of lO^^G. 

In the case of a homogeneous magnetic field a more direct single frequency measurement 
is possible. The basis of this technique is to exploit the parity and symmetry properties 
that such a field induces in the CMB. Given a measurement of the Stokes parameters, as 
a function of position on the sky, this data can be decomposed into a sum of generalized 
spherical harmonics. Thus 



^(n) = ^a^F;„(n); 

Im 

{Q±iu){iv) = Y.4^±2yu^) (1) 



Im 



where T denotes temperature and the spin ±2 spherical harmonics, ±2^/™,, are used to 
preserve the behavior of Q and / under coordinate rotations. Here we follow the notation 
discussed in |]13|; for an alternative, but equivalent approach see fl^. Circular polarization 



can be decomposed as a scalar, but is not expected to occur in the CMB [ p^ . 
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Under parity inversion 
This motivates us to define 



-lY_gYem and thus 2yem ±-2 are parity eigenstates. 
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lYaf^ and a£, 



and 



lYaf^ under parity inversion. In an isotropic universe, 
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SO that af^ — > 

cross correlations between the B polarization and the temperature and E polarizations are 
forbidden as this would imply noninvariance under parity. A homogeneous magnetic field, 
however, is maximally parity violating and therefore one would expect one of its primary 
signals to be such a cross correlation. Motivated by this expectation, we consider in some 
detail the effect of a homogeneous magnetic field on CMB polarization. 

Polarization of the CMB is generated in the optically thin last scattering surface by 
quadropole fluctuations in temperature. Scalar fluctuations result in polarization that is 
purely of the E sort, while tensor fluctuations result in B and E polarization. The presence 
of a magnetic field, however, will create B polarization that is correlated with the measured 
temperature spectra. This happens as the the B and E polarizations are Faraday rotated into 
each other. Following |[T^, we can estimate the extent of this coupling by taking the optical 
depth of the last scattering surface to be approximately 1. For a homogeneous magnetic 
field we find 



2.3'= 



( 



10-9G/ V30GHz 



-2 



cos^^ 
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where Q is the angle of the magnetic field with respect to the direction of propagation. This 
effect can be used to determine the magnetic field by comparing observations at different 
frequencies |T^. As this B-polarization signal is expected to be small, however, one would 
like to minimize experimental noise as well as systematic uncertainties by comparing it to 
quantities at the same frequency and with greater amplitude. Given our discussion above 
{(ie*m'(^fm) represents just such a comparison. We would therefore like to examine this 
quantity in greater detail. 

The evolution of radiation in a perturbed Friedman-Robertson- Walker universe can be 
studied in the linear regime. In this case each plane wave perturbation can be considered 
separately, with an average over the k's being done at the end of the calculation. Defining 
/i = n ■ k where n is the line of sight we can write down the Boltzmann equation in the 



synchronous gauge as follows 



At + ikfiAx 



Aq + ik^Aq 
Au + ikfiAu 

n 



+ k 



-At + Ato + ifiVb + -P2(/u)n 



-AQ + -[l-Q2{fi)]U 



- kAu - 2ujbAq 

At2 + Aq2 + Ago 



2ujbAi 



(4) 

(5) 
(6) 



where At-, Aq, Au are the brightness functions for T, Q and U, which can be expanded 
in multipole moments defined such that A(r, k,^) = I]£(2^+ 1)(— i)^A^(r, k)Pf{n) {Pe{n) is 
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the Legendre polynomial of order i). Derivatives are taken with respect to the conformal 
time r, k is the differential optical depth for Thomson scattering k = {a/ aQ)nf.aT and h, rj, 
Vb are sources of temperature anisotropics from metric perturbations and baryon velocities. 
We shall restrict ourselves to scalar perturbations. The effect of Faraday rotation comes in 
through the coupling between polarization components, wb = k^ujQCOs9 where 



ujq = .06 



Bo 
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we see that Aq (A(/) give 



and we have assumed that Bq lies along the z axis. From Eq. 
us the amplitude of the E (B) component of the polarization. 

Before we proceed, we shall now discuss two important approximations. Firstly we will 
consider an isotropic universe supporting a homogeneous magnetic field. As argued in 
this is not entirely correct, such a field can only be supported by the existence of a globally 
anisotropic term. However, as we shall see, most of our statistic will rely on the large £ 
behavior of the cross correlation spectrum and will therefore be insensitive to the large scale 
behavior of the anisotropy. Clearly a more detailed calculation should include such terms 
(for a good description see fl^) but for the purpose of this letter we shall not pursue it. 
Secondly, we are interested in studying Bq in the regime where it is competitive with other 
constraints from the CMB, notably ||^, and consequently ub is sufficiently small that we can 
drop it from Eq. This means that the only source term for the correlated B polarization 
is through a rotation of E polarized light by ~ 2°; it is expected to be small compared to 
E and the coupling of B into E even smaller. Thus the effect can be applied only to the 
evolution of B polarized light, and there only as a scale factor of Bq cos 6' in the source term. 

This represents a great simplification in the calculation. As the angle of the magnetic 
field with respect to the observer appears only as a scale factor, we can sidestep the issue of 
the relative orientation of the wave vector k and Bq, using it only in the angular averages 
taken at the end of the calculation. With this simplification we find the cross correlation 
due to scalar perturbations to be 
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where 
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P{k) is the initial power spectrum of scalar perturbations, x = kijQ — r), a = [h — 6fj)/2k'^, 
je{x) is the spherical Bessel function of order i, and g{T) = ke'^. Physically, g{T) is a 
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visibility function, whose peak defines the epoch of recombination, which gives the dominant 
contribution to observed anisotropics. Note the double (?(r) integral in the polarization term; 
this appears as during last scattering the T anisotropics are the source of E polarization 
which in turn is the source of B polarization. 

By expressing cos 6 in terms of Yi^ and applying the three Y^^m formula, Eq. |^ can be 
further simplified as 



4£2 - 1 



1/2 



4(£+ 1)2-1 



1/2 
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The form of our result illustrates the SO (2) symmetry of the model which we are considering 
||r8| . To see why this is so, consider the general case of a^m's that are generated by a model 
that is cylindrically symmetric about the z axis. Then the symmetry transformations are 
R{(p)o,em = G^'^'^aem and so {ae'm'dem) = ^mm'C^^^ as terms with m ^ m' are not invariant 
under cf) rotations. Furthermore, as we have found a nonzero correlation between coeffi- 
cients of opposite parities, this indicates that the axis of symmetry has an orientation. Our 
magnetic field lies along a definite axis and it points a definite way along it. The effect of 
a Bq pointing in the x or y directions can be similarly obtained by replacing the cos^ in 
Eq. H with X = sin 9 cos (p and y = sm9 sin 0, yielding equations for Cj^i ^ in the x and y 
directions. 

Thus we have discovered measurable quantities that directly scale with the magnitude 
and direction of the magnetic field. Experimentally, ^(xi) can be determined by ar- 
bitrarily choosing a basis on the sky and then using appropriate combinations of the a^'s 
and to form estimators in the x, y and z directions. We have calculated ; using the 



CMBFAST code of |]T9|. We shall restrict ourselves to a standard cold dark matter (CDM) 
universe to quantify the effect of the magnetic field; the results appear in Figure |^ and where 
we normalize by the scale of the magnetic field and plot //(^o- For comparison, we also 
plot Cjf. Note that if the last scattering surface were taken to be infinitely thin, then 
g^r) = 6{t — r*) and C^^/wq and Cff would be identical up to the difference introduced 
by integrating over two bessel functions whose orders differ by one. As polarization is gen- 
erated in the optically thin last scattering surface, the visibility function is quite sensitive to 
reionization. Thus we have also calcuated ^ for a CDM universe with an optical depth 
of K = 1 and we present the results in the bottom panel of Figure |1[ Note that reionization 
will supress power in the various cross correlations. However the broader support of g^r) 
leads to a slight increase in amplitude of Juo relative to Cjf. 

Let us now turn to the observability of such a signal. For the purpose of this paper we 
shall consider a full sky survey such as would be performed by a satellite. To estimate the 
error in determining C'^^ we follow the analysis given in ||14[. Being careful to include the 
m coeficients in Eq. [1^ we find 
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(11) 



where a is the full width at half maximum of the beam, and Wp are measures of the 
experimental sensitivity, independent of pixel size and we have assumed full sky coverage. 
We have also approximated 
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FIG. 1. The top panel shows a comparison of i{£ + l)Cjf /{2'k) (soUd line), 

t{l + l)Cj4//(2vru;o) (dotted line) and t{i + l)C|^^^^/(27rwo) (dashed line) in a standard CDM 
cosmology, while the bottom panel shows such a comparison in a reionized universe with k = 1. 
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which is vahd for large ts. Due to the high cosmic variance and low signal at small i we 
will want to base our measurement of the magnetic field on correlations with i > 50 where 
this approximation is justified. 

Following 1^ we can define a function for the goodness of fit of a theory given an 
observation: 



X^(^o) = E 



r^TBexp ^TBth(, , i^TBexp riTBtht, , Xl 



(13) 



where the sum is taken from £ > 50 and we assume that we know all other cosmological 
parameters. We hope the true value of the magnetic field will minimize the its 
sensitivity to Bq will tell us how well we can constrain it. We are interested in determining 
the minimum value of Bq to which we are sensitive. In this limit we expect noise to dominate 
over C^^ and C^^ as they scale as Bq and Bq respectively, so we can neglect these terms in 
Eq. 11. 

Following pil 



an approximate l—a error is given by 
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In Figure y we show two contour plots for the expected precision in {]1^^){^q;^) 
for a range of experimental parameters. Here we take wp = 2wt WM- For low frequency 



detectors (z/q < 50Ghz), the top plot shows that we can get comparable limits to that of 
for wt < 10^^^. Such frequencies are currently only accessible through HEMT detectors 
where the sensitivity has yet to reach such high levels. With low frequencies there is the 
additional problem of having to consider large beamwidths. For example an optimistic, 
space based 30GHz receiver can achieve at most 15' resolution. The planned MAP mission 
3g will have a 30GHz with a = 42' and w^^ ~ 4.7 x 10"^^ the Planck mission |2§] will have 
a 31GHz detector with a = 30' and slightly smaller wt- This means we will get a constraint 
of (j{Bo) — few X 10~*Gauss from both of them. 

For high frequency measurements one is faced with the term, i.e. the higher the 
frequency the smaller the effect. In the bottom plot of Fig ^ we have consider a lower 
range of wt appropriate for these measurements, and labeled isocontours in steps of 0.1. 
The important thing to note is that current bolometer technology is reaching these levels 
of sensitivity. If we consider the example of the 150GHz detector on the proposed Planck 
mission, we find that it may be possible to get a constraint of a{Bo) — 10~^, an improvement 
over the limits set in [Q. With the rapid advances in both high frequency HEMT and 
bolometer technology it is conceivable that one may do even better (for a realistic assessment 



of future prospects see ||2^ ) 
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In this letter we have presented a new method for constraining a homogeneous, pri- 
mordial, magnetic field. In describing the technique and its potential sensitivity, we have 
restricted ourselves to a full sky measurement of the CMB anisotropy and polarization. 
Clearly this is a simplification. A full analysis should include different experiments with 



varying degrees of sky coverage We have also focused on scalar perturbations in a 

CDM universe. The inclusion of gravity waves and vorticity will generate a non-negligible 
C^^ which may increase the cosmic variance of our statistic and necessarily reduce its sen- 
sitivity. Naturally the inclusion of global anisotropy, the effects of foreground sources and 
different thermal histories may also change our results. For example, for the reionized uni- 
verse considered in Figure |I] the constraint is worse by a factor of two as compared to the 
sCDM case. Again a full analysis should include a wide range of cosmological parameters 
and scenarios. 

This method has two appealing features which we restate. Firstly all current scenarios 
of structure formation assume statistical isotropy which necessarily leads to a zero cross 
correlation between B type polarization and temperature anisotropy. It is only through the 
existence of magnetic field (or some form of global anisotropy) that such a cross correlation 
can be induced. Secondly the constraint derived in P] was strongly dominated by the large 
angle anisotropics. The authors pointed out that a good measurement of the temperature 
anisotropics already exists on these scales and future measurements will do little to improve 
it. This is not the case of our method. Clearly, the better the measurements become the 
better the constraint of Bq will be. 
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